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Nonclassicality cannot be a single-observable property since the statistics of any quantum ob¬ 
servable is compatible with classical physics. We develop a general procedure to reveal nonclassical 
behavior from the joint measurement of multiple observables. This requires enlarged spaces and 
couplings with auxiliary degrees of freedom so it must be followed by an inversion procedure to infer 
system properties from the observed statistics. In particular this discloses nonclassical properties of 
standard examples of classical-like behavior, such as SU(2) and Glauber coherent states. Moreover, 
when combined with other criteria it would imply that every quantum state would be nonclassical. 
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I. INTRODUCTION 

Nonclassicality cannot be a single-observable property 
since within classical physics it is always possible to re¬ 
produce exactly the statistics of any quantum observable. 
For example, regarding optics, a classical wave with def¬ 
inite complex amplitude is infinitely squeezed and max¬ 
imally subPoissionian. Moreover, we can produce classi¬ 
cal waves whose intensity replicates any discrete photon 
number distribution as accurately as desired. Thus, non¬ 
classical effects must be found in the joint statistics of 
multiple observables, specially if they are incompatible 

DP- 

Joint measurements require the coupling of the system 
with auxiliary degrees of freedom. This means that in 
general the measurement must be followed by some kind 
of data analysis or inversion procedures to extract the in¬ 
formation of the system from the observed statistics [2]. 
We will show that when this is done in classical physics 
the result of the inversion is always a true probability dis¬ 
tribution of the system variables. However, in quantum 
physics this is not always the case and the result of the 
inversion can be incompatible with classical statistics. 

This provides a very general program that includes 
as particular cases classic demonstrations of nonclassi¬ 
cal behavior, such as coincidence detection and quantum 
tomography. The retrieved system information is usu¬ 
ally expressed in terms of the Wigner W or Glauber- 
Sudarshan V functions, that actually intend to represent 
the joint distribution of two complementary field quadra¬ 
tures [31 03 . 

It can be expected that this program will provide 
meaningful results when gathering information about 
complementary observables. But it can be as well ap¬ 
plied successfully when addressing the redundant mea¬ 
surement of the same observable, such as photon num¬ 
ber. Besides providing a new perspective on the subject, 
this approach discloses nonclassical properties for states 
with classical-like V distribution, that otherwise are uni¬ 
versally considered as classical light, with few exceptions 
EH23- Moreover, when combined with other nonclassical 
criteria it would imply that every quantum state would 


be nonclassical. 

The details of the program are presented in Sec. II, in¬ 
cluding that in classical physics the result of the inversion 
is always a true probability distribution. The program is 
applied to a qubit via the measurement of complemen¬ 
tary observables in Sec. III, showing nonclassical spin 
features for SU(2) coherent states. In Sec. IV we go be¬ 
yond complementarity considering the redundant mea¬ 
surement of the photon number. Finally, in Sec. V we 
show a nonclassical effects displayed by Glauber coherent 
states after the joint measurement of phase and number. 


II. PROGRAM 

We consider the simultaneous measurement of two 
compatible observables A' and Y which is intended to 
provide information about two system observables X and 
Y [8j. The statistics of the measurement px,y{x,y) is 
arranged so that each variable x , y refers to the corre¬ 
sponding observables A, Y with marginal distributions 

px (x) = X p x x ( x >y)’ py ( v ) = X p x ’ y y) ■ 

y x 

(!) 

We have assumed a discrete range for x and y without 
loss of generality. The joint measurement usually requires 
the coupling of the system with auxiliary degrees of free¬ 
dom that must be followed by some kind of data analysis 
in order to extract information about X and Y in the 
system state. To this end the only hypotheses we make 
is that it is possible to infer the exact distributions pz{z ) 
from their observed counterparts pz{z), with Z = X,Y, 
in the form 

Pz{z) = YlP / z{z,z')pz{z'), p z = pzPz, (2) 

z' 

where \iz is a matrix with matrix elements pz{z, z') com¬ 
pletely known as far as we know the measurement being 
performed. These are the non ideal invertible measure¬ 
ments introduced in Ref. [2|. 
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The key idea is to extend this inversion from the 
marginals to the complete joint distribution in the form 

Px,v(x,y) = ^2 Px{x,x')p, Y (y,y')p x , Y (x' ,y'), (3) 

x',y' 

or in matrix form 

PX,Y = PxPx,yPy, (4) 

where the superscript t denotes matrix transposition. 

The retrieved distribution px, Y (x,y) can be referred 
to as Wigner measure since by construction px, Y {x,y) 
provides the correct marginals for X and Y as a typical 
property of Wigner functions. In particular, when A' and 
Y are field quadratures the result of the inversion is the 
classic Wigner function [2]. 

Classical domain. - Next let us show that in the classi¬ 
cal domain this procedure always leads to a bona fide dis¬ 
tribution px,Y( x ,y). Classically the state of the system 
can be completely described by a probability distribu¬ 
tion P(A), where A is a point in the corresponding phase 
space. We may say that A represent the ontic states, 
while P( A) are the epistemic states :9]. The observed 
joint statistics is always of the form 

Px,y{x,v) = J dXP(X)X(x\X)Y(y\X), (5) 

where Z(z\\) for Z = X,Y are the conditional probabil¬ 
ities that the observable Z takes the value 0 when the 
system state is A. Imposing that Eq. (|2j) holds 

Z(z\X) =^2pz(z,z')Z(z'\X), (6) 

z' 

we readily get from Eq. ([ 3 ]) the actual joint distribution 
for X and Y 

Px,y( x , y) = J d\P(X)X(x\X)Y(y\X). (7) 

This is to say that after the form © the inversion works 
equally well for the joint distribution as for the marginals. 
Therefore, the lack of positivity or any other pathology of 
the retrieved joint distribution px. Y (x,y) is a signature 
of nonclassical behavior. 


III. QUBIT COMPLEMENTARITY 

The simplest example holds in a two-dimensional 
Hilbert space considering two complementary observ¬ 
ables represented by the Pauli matrices X = a x and 
Y = cry, with outcomes x : y = ± 1 . The most general 
system state is represented by the density matrix in the 
system space H s 

P = \ ( CT o + s ■ er ), (8) 


where <7o is the identity, er are the three Pauli matrices, 
and s = tr(per) is a three-dimensional real vector with 
s < 1. The exact statistics are 

Px(x) = \ (1 + xs x ) , p Y (y) = ^ (l + ys v ). (9) 

For example we consider that Y is measured directly 
on the Jis variables by projection on its eigenvectors 
|y = ±1). The simultaneous measurement of X can be 
achieved by coupling the system space T-L s with auxiliary 
degrees of freedom via the following unitary operator 

uni 

U = V + \x = l)(x = 1| + V-\x = —Y){x = — 11, (10) 

where V± are unitary operators acting solely on and 
\x = ±1) are the eigenvectors of X = a x . Denoting 
by | a) the initial state in the transformed states are 
|a±) = V±\a) G ' H m ■ They are not orthogonal in general, 
and their overlap (a+|a_) = sin<^> is assumed to be a 
positive real number without loss of generality. The most 
simple measurement in TL m is described by projection on 
the orthonormal vectors |i) 

|i = 1> = !«+>-»■> 11«->), 

|i = _1> = ^(“ sil, 5 |o + > + c “5 |a - ) )' (11) 

with (l|a+) = (—l|a_) = cos(^>/2), (l|a_) = (—l|a+) = 
sin(</>/2). 

With all this, the joint statistics for the simultaneous 
measurement is 

Px, Y (x,y) = ^ (1 + xs x cos (j) + ys y sin </>), (12) 

with marginals 

px(x) = ~ (1 xs x cos (jP ), p Y {y) = \ (1 + ys y sin(f) , 

(13) 

which can be easily inverted in the form 

Pz = PzPz , p z {z,Y) = \ (1+ — ) , (14) 

2 V Pz J 

with r/x = cos <fr, rjy = sin<^, and z. z' = ±1. Fi¬ 
nally, the inversion applied to the joint statistics px,Y = 
PxPx.vPy leads to 

Px,v{x,y) =-{l+xs x +ys y ). (15) 

This is not a true Wigner function in the sense that the 
information about the a z observable is absent. 

Clearly Px,y(x,p) can take negative values if |s| > 
l/-\/2- The states presenting the largest negativity are 
pure states with |s x | = |s. y | = l/y/2 and s z = 0. The sys¬ 
tem space for the qubit is the unit Bloch sphere. Relying 
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on the spherical symmetry we may say that any state 
with |s| > l/v 2 is actually nonclassical after a suitable 
choice of the X , Y observables. This is that the 65 % of 
the volume of the Bloch sphere is occupied by nonclas¬ 
sical states. This result is quite relevant since it is often 
argued that all states in a two-dimensional space are clas¬ 
sical regarding spin variables, as far as all them have a 
well-behaved Glauber-Sudarshan SU( 2 ) V function jTlj. 
Actually, all the pure states are SU( 2 ) coherent states 
and have the maximum negativity. 


IV. NUMBER-NUMBER MEASUREMENTS 

Next the system is an electromagnetic field mode. This 
is mixed with vacuum at a lossless 50 % beam splitter. 
The observables to be measured X , Y are the number 
of photons ni and 772 registered at the two output ports. 
Both are intended to provide information about the num¬ 
ber of photons in the system state p, this is X = Y = n. 
The observed joint statistics is (for simplicity we drop the 
subscripts in p and p since there is no risk of confusion) 

P(ni,n 2 ) = ^^( ni + i n 2 ^p(n 1 +n 2 ), ( 16 ) 



FIG. 1 : Chessboard with the p(ni,ri2) values for a number 
state |n = 7 ). The blue squares represent negative p(ni,ri2) 
values, the darker the more negative. 


From the above relations we get that if V(a ) is a classical- 
like distribution so will be £>(771,772), provided it exists. 


with identical marginals for both variables 77 = 77 i, 772 

P( n ) = fl^k (^W), (17) 

k—n V / 

where p[n) = ( 77 |p|n) is the true number statistics of 
the system state. The inversion of the marginals can be 
carried out via a transformation readily similar to Eq. 

p(n) = (-2)»f;(M(-l) fc p(fc). (18) 

k=n ' ' 

When the inversion is applied to the two variables in the 
joint distribution p(rii, 772) via Eq. © we get a distribu¬ 
tion p(m, 772) that after a little algebra can be expressed 
as 


p(ni,n 2 ) = Y ^ Ul + , n2 ,t| fc ^ (~ 1 ) fc ^( n i+ n2+fc )> ( 19 ) 

^0 ni]n2lkl 


with a clear symmetry between 771 and 772. This adopts a 
quite simple expression if we use the Glauber-Sudarshan 
distribution V(a) in the form 

£>( 771 , 772 ) = — p —7 f d 2 aP(a)\a\ 2ni+2n2 e~ 2 \ a \ . (20) 

77 i! 77 2 ! J 


In particular for m = 772 = 0 we get that £>( 0 , 0 ) is the 
Wigner function W (a) at the origin a = 0 as propor¬ 
tional to the mean value of the parity operator 


p(0,0) = |VF(a = 0). 


Number states. - A readily example is provided by the 
photon-number states I77). For example for the one- 
photon state 77 = 1 we have £>( 0 , 0 ) = — 1 , p(l, 0 ) = 
£>(0,1) = 1, while all the other vanish p(ni 7^ 0,77-2 7^ 
0 ) = 0 . For states with larger number of photons we 
get results as illustrated by the chessboard in Fig. 1 for 
1 77 = 7 ). The blue squares represent negative values for 
p(n 1,772), the darker the more negative, with an absolute 
minimum at 77 i = 712 = 2 with £>( 2 , 2 ) = — 210 . Actu¬ 
ally, the absolute value of the most negative £>(771,772) 
increases heavily with the number of photons 77, with 
rather strong positive-negative oscillations. 


Photon added thermal states. - As a more realistic ver¬ 
sion of the number states we may consider the photon- 
added thermal states p oc a^ k p t a k P 2 |, where p t is a 
thermal state, with photon-number distribution 

^ = m 

and £>(77 < k) = 0 , where n is the mean number of pho¬ 
tons of the thermal state. It can be easily seen that for 
the single-photon-added thermal state k = 1 the only 
negative value £>(771,772) < 0 holds for 771 = 772 = 0 with 
p( 0 , 0 ) = — l/( 2 fi + l) 2 , as illustrated in Fig. 2 . 


Squeezed vacuum. - As a further example we may con¬ 
sider the squeezed vacuum state with photon-number dis¬ 
tribution that contains only contributions for even num¬ 
bers [ 13 ] 

77! tanh” r 
^ n 2 n [n/2)\ 2 coshr 


( 21 ) 


(23) 
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FIG. 2: Chessboard with the p{n\,ri 2 ) values for a single¬ 
photon-added thermal state k = n = 1. The blue squares 
represent negative p{n\,ri 2 ) values, the darker the more neg¬ 
ative. 



FIG. 3: Chessboard with the p(ni,ri 2 ) values for a squeezed 
vacuum state with n ~ 0.1. The blue squares represent neg¬ 
ative p(n i, 722 ) values, the darker the more negative . 


with p(n) = 0 for odd n, where r is the squeezing pa¬ 
rameter and n = sinh 2 r is the mean number of photons. 
The most negative result p(ni,ri 2 ) < 0 holds always for 
p( 0,1) = p( 1,0) = —fi as illustrated in Fig. 3 for r = 0.3 
and n ~ 0.1. Numerically we have shown that the in¬ 
version for the joint distribution works only for n < 0 . 1 , 
while for the marginals it still works about this value. 


V. PHASE-NUMBER MEASUREMENTS 

Let us address again a measurement scheme with two 
photon-number detectors placed at the output ports of 


a 50 % lossless beam splitter providing a joint statistics 
p' ( 211 , 712 ). However, now the beam splitter mixes the 
system state p with a Glauber coherent state |/3). In 
this case, the variables of interest are not the number 
of photons recorded in each detector 711 , 2 , but the total 
number N = n\ + 71-2 and the normalized number differ¬ 
ence m = (n 1 — 1 ^ 2 )/(ni + 712 ) with m £ [—1,1]. The 
measured joint distribution for these variables is 

( "1 ff) 1 — 7?7 \ 

n\ = N —^—, n 2 = N —-— J , (24) 


where N = 0,1,. ..,00, while the range for m is the union 
of all the ranges allowed for each N, this is m = —1, — 1 + 
2/N, ..., 1, and we consider just m = 0 for N = 0. 

The marginal statistics for the total number p(N) can 
be regarded as providing information about the photon- 
number distribution p[n) = (n\p\n) of the system state 
p. On the other hand, the marginal distribution p(m) is 
a suitable version of the cosine of the relative phase. For 
the sake of simplicity we consider that p{m) requires no 
further data transformation. This is consistent with the 
understanding that there is no operator for the single¬ 
mode phase, and that this is a somewhat fuzzy variable 
represented by positive operator measures [HI ITS] . In 
particular, the variable m bears some resemblance with 
relative-phase approaches introduced via homodyne de¬ 
tection [IMT5] . 

The marginal for the photon-number sum is 

N -N—n 

p( N ) = p( N ’ m ) = e ~ n ( 25 ) 

me# n=0' ' 


where n=\j 
herent state 
in the form 


| 2 is the mean number of photons of the co- 
0 ). The above relation (25) can be inverted 


77 / — \ti _ A / 

N—0 V J 


(26) 


which is essentially the same transformation (251 after 
replacing n by — n. The program ends applying the in¬ 
version transformation to the complete joint distribution 


p(n,m) = e ^ j— 


l-N 


N—0 


(n - N)\ 


■p(N , m). 


(27) 


A readily demonstration that this can disclose non- 
classical behavior holds when considering the particular 
values n = 1 and m = 0 , taking into account that m = 0 
is not included for N = 1 


p(l, 0 ) = —ne n p( 0 , 0 ) = —ne n p'{ 0 , 0 ) = — np(0). (28) 

Therefore, whenever p(0) 7 ^ 0 and fi ^ 0 we have 
p(l, 0 )< 0 . 

This result may be ascribed to same kind of effective 
entanglement between the N and m variables since the 
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FIG. 4: Chessboard with the p(N, m ) values for p in the vac¬ 
uum state mixed in a 50 % beam splitter with a reference 
coherent state with n — 1. The blue squares represent nega¬ 
tive p(N, m) values, the darker the more negative. 


The result of the inversion (27) is illustrated in Fig. 4 for 
n = 1 where the most negative value is P(N = 1 , m = 
0 ) = - 1 . 

The fact the coherent states display nonclassical be¬ 
havior for number-phase variables was already noticed in 
Refs. 5j Ifi. 


Every quantum state is nonclassical . - To some extent 
the result in Eq. (28) is the dual of the nonclassical Lee 
criterion in Refs. 19], where maximal non classicality is 
obtained by the lack of vacuum contribution. In other 
words, the union of both criteria would imply that all 
quantum states are nonclassical. 


allowed m may be different for different N. This is a 
purely quantum-phase effect since in the classical case 
there is no relation whatsoever between the intensity and 
the relative phase that can take any value between -1 
and 1 for any value of the intensity. This is reflected 
by most approaches to the quantum relative phase [El- 
Equivalent results are obtained for other definitions of 
the m variable such as m = n\ — n 2 or even m = n\. 


Nonclassicality of Glauber coherent states. - In partic¬ 
ular Eq. (28) discloses non classical properties for the 
Glauber coherent states |a), which are widely assumed 
as classical-like states. The most simple example is the 
vacuum state \a = 0) with 


VI. CONCLUSIONS 


Summarizing, we have addressed an universal proto¬ 
col to disclose nonclassical behavior via joint measure¬ 
ments of more than one observable. This includes pre¬ 
vious criteria and admits many other possibilities under 
one and the same framework. In particular this reveals 
nonclassical properties for standard examples of classical- 
like behavior, such as SU(2) and Glauber coherent states. 
Moreover, when combined with other criteria implies that 
every quantum state is nonclassical. 


and 


p' (ni,n 2 ) 


(»/2)" 1+ " 1 2 3 4 5 c _ fi 

ni!n 2 ! 


p(N , to ) 


(n/2) N 

N (1+m) | N(l—m) | 

2 ‘ 2 ’ 
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